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ABSTRACT 

We present a spectral line survey of the C-rich envelope CIT 6 in the A 2 mm 
and 1.3 mm bands carried out with the Arizona Radio Observatory (ARO) 12 m 
telescope and the Heinrich Hertz Submillimeter Telescope (SMT). The observa- 
tions cover the frequency ranges of 131-160 GHz, 219-244 GHz, and 252-268 GHz 
with typical sensitivity limit of Tr < 10 mK. A total of 74 individual emission 
features are detected, of which 69 are identified to arise from 21 molecular species 
and isotopologues, with 5 faint lines remaining unidentified. Two new molecules 
(C4H and CH3CN) and seven new isotopologues (C^^O, ^sSiCz, ^^SiO, =^°SiO, 
^^CS, C'^^S, and C'^'^S) are detected in this object for the first time. The column 
densities, excitation temperatures, and fractional abundances of the detected 
molecules are determined using rotation diagram analysis. Comparison of the 
spectra of CIT 6 to that of IRC+10216 suggests that the spectral properties of 
CIT 6 are generally consistent with those of IRC+10216. For most of the molec- 
ular species, the intensity ratios of the lines detected in the two objects are in 
good agreement with each other. Nevertheless, there is evidence suggesting en- 
hanced emission from CN and HC3N and depleted emission from HCN, SiS, and 
C4H in CIT 6. Based on their far-IR spectra, we find that CIT 6 probably has a 
lower dust-to-molecular gas ratio than IRC+10216. To investigate the chemical 
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evolution of evolved stars, we compare the molecular abundances in the AGB 
envelopes CIT6 and IRC+10216 and those in the bright proto-planetary nebula 
CRL618. The implication on the circumstellar chemistry is discussed. 

Subject headings: ISM: molecules — radio lines: stars — line: identification — 
stars: AGB and post-AGB — stars: circumstellar matter — stars: individual 
(GIT 6) — surveys 



Introduction 



The late stages of stellar evolution from the asymptotic giant branch (AGB) to planetary 
nebulae (PN) are now recognized as an active period of chemical synthesis of molecules. The 
detection and analysis of millimeter wave molecular emission lines are fundamental to the 
understanding of the physical conditions and chemical processes leading to chemical synthe- 
sis. Due to the rapid evolution of the star, the changing physical conditions, including dust, 
stellar winds, shock waves, UV emission and X-rays from the central star, etc. play different 
roles in circumstellar chemistry. This leads to corresponding different circumstellar chemical 
compositions in different evolutionary stages. The envelopes around G-rich stars, with their 
enhanced carbon abundance, provide a perfect cradle for molecule f ormation. Hitherto, more 
than 60 molecular species have been detected i n G-star envelopes (lGlassgoldlll996l : lOlofsson 
19971 : iGernicharo et al.ll200Cll : IZiurys et al.ll2007l ) , most of which were discovered through their 
rotational lines at millimeter wavelengths. 

Recent improvement in telescope design and receiver performance enable us to detect 
new molecular emission with a higher sensitivity, and thus with the possibility of shedding 
new light on circumstellar chemistry. The most frequently investigated G-star envelope is 
IRG+ 10216, which is one of the richest molecul ar sources in the sky. Several molecular line 
surveys have been presented for this object (see IGernicharo et al.ll2000l : iHe et al.ll2008l . and 
the references therein), which was found to harbor extremely abundant carbon chain and 
metal-containing molecules. IRG+10216 has been frequently used as a standard reference 
for the chemical compositions of late-type stars. This inevitably invites the issue whether 
IRG+10216 is a chemically unique late-type star. To settle this question, we require system- 
atic surveys of molecular line emission from other G-star envelopes. In the present study, we 
report a spectral-line survey of the G-star envelope GIT 6 at millimeter wavelengths. This 
allows us to compare the similarity and difference in the chemical compositions between the 
two G-star envelopes. 



GIT 6 (RWLMi, GL 1403, IRG+30219, IRAS 10131+3049) was first discovered during 
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the Caltech 2- iim sky survey an d was among the 14 very red infrared-bright optical-faint 
sources found (jUhich et al.lll966l ). CIT6 is characterized by its very low color temperature, 
implying that the star is surrounded by a very thick dust envelope and has been identified 
as a long-period var iable with a period of ab out 628 days (jAlksnislll995l ). From the period- 
luminosity relation, ICohen fc Hitchoru (119961 ) estimated the distance of CIT 6 to be 400 ± 
50 pc, ^ yhich is slightly more dis tant than IRC-l -10216, which has a d istance of between ~ 
120 pc (iGroenewegen et al.lll998l ) and ~ 150 pc (ILucas fc Guelinlll999l ). CIT 6 is believed to 



bemore evolved and has a lower mass loss rate compared to IRC-l-10216 (see iFukasaku et al. 
19941 . e.g.). 



The large polarization found in the visible and infrared wavelength s implies that the 



distribution of c ircumstellar material around the star is asymmetric (iKruszewskil Il968 



Dyck et al.l Il97ll ). Multi- wavelength imaging observations have been performed to study 
the structure of the nebula around CIT 6. The optical images obtained by the Hubble Space 
Telescope (HST) and the near-IR images of CIT 6 obtained by the Keck-I telescope have 
revealed a bipolar du st envelope and an elongated component with time-variable asymmetry 
f Monnier et al.l 2000). Several sca ttering arcs were revealed by the HST-NICMOS imaging 
polarimetry (jSchmidt et al.ll2002l ). These arcs are ne arly concentric and extend to large 
stellar radii. Mid-IR images of CIT 6 were obtained by iLagadec et al.l (120051 ) using the ESO 
3.6-m telescope. A cometary-like feature was revealed in their 9.7/im image. 

There have been several observations of molecular lines in CIT 6 at millimeter wave- 



length. iHenkel et al.l (119851 ) reported observations of a few molecular lines in CIT 6 and 
IRC+10216 between 18 and 150 GHz. They found that relative abundances of observed 
molecules in the two sources hav e no significant differences. Using the Nobeyama 45 m radio 
telescope, iFukasaku et al.l (119941 ) observed a few transitions in the frequency ranges between 
39-47 GHz and 85-91GHz in a sample of evolved stars including CIT 6 an d found that the 



abundance of HNC increases with the evolutionary stage of the stars. iBujarrabal et al. 



(119941 ) presented observations of 10 molecular transitions in C-rich and 0-rich circumstellar 
envelopes including CIT 6 with the IRAM 30 m radio telesco pe at 1.3 mm, 2 mm , and 3 mm 
windows. A recent molecular line survey was presented by IWoods et al.l (120031 ) using the 
SEST 15 m and Onsala 20 m telescopes. They found that CIT 6 stands out from the other 
C-rich envelopes due to its high CN/HCN ratio and low HNC/HCN ratio. To date, the 
molecular species positively detected in CIT 6 at millimeter wavelengths include CO, ^^CO, 
CN, i^CN, CS, SiO, SiS, 29Si32S, CaH, SiCa, HCN, H^^CN, HNC, C3N, HC3N, HC^^CCN, 
HCC^^cN, and HC5N. 



In this paper, we present the first systematical line survey of CIT 6 at the 2 mm and 
1.3 mm windows, using the Arizona Radio Observatory (ARO) 12 m telescope and the Hein- 
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rich Hertz Submillimeter Telescope (SMT). The observations are described in Sec. 2. In 
Sect. 3 we present the identifications and abundance calculations of the detected molecular 
species. In Sect. 4 we discuss the implication of our findings on circumstellar chemistry. The 
conclusions are given in Sect. 5. 

2. Observations and data reduction 

The observations were carried out in beam switching mode with an azimuth beam 
throw of 2' during the period from 2005 April to 2006 September. Pointing was checked 
by observations of a planet every 2 hours. The 131-160 GHz (2 mm window) spectra were 
obtained with the ARO 12 m telescope at Kitt Peak, and the 219-244 GHz and 252-268 GHz 
(1.3 mm window) spectra with the SMT 10 m telescopes on Mount Graham, Arizona. The 
2- and 1.3- mm dual-channel SIS receivers were employed, operated in single sideband dual 
polarization mode with a typical image rejection ratio of > 18 dB. At the ARO 12 m, the 
spectrometer back-ends were two 256-channel filter banks (FBs) with a channel width of 
1 MHz and a millimeter autocorrelator (MAC) with 3072 channels and 195 kHz per channel. 
The spectrometers utilized at SMT were a 2048-channel acousto-optical spectrometer (AOS) 
with a channel width of 500 kHz and 1024-channel Forbes Filterbanks (FFBs) with a channel 
width of IMHz. The system noise temperatures were typically 150-400 K at 2 mm and 400- 
700 K at 1.3 mm. The temperature scales at the ARO 12 m and the SMT, and T^, were 
obtained using standard vane calibration. A ~ 15% calibration error was obtained from a 
few strong lines detected in different spectrometers and different epochs. The main beam 
brightness temperatures were derived through Tr = T^/rj^^ and Tr = T^/rjmb, for the 12 m 
and the SMT data, respectively, where r]^ is the corrected beam efficiency (~ 0.75; see the 
ARO 12m manual for its definition) and rj^b the beam efficiency (~ 0.7). Both telescopes 
have a large beam size, covering the whole emission region. 

The CLASS software package in GILDAS was used to reduce the spectra. After 
discarding the bad scans which are seriously affected by bandpass irregularities, we co-added 
the spectra from individual scans. The baseline was fitted with a low-order polynomial. In 
order to improve the signal-to-noise ratios, the spectra were smoothed and rebinned by a 
factor of 3, yielding a typical rms noise temperature of < 10 mK in main beam brightness 
temperature units and a spectral resolution of ~ 4km/s. Since each line has been observed 
by two different spectrometers, we can reduce the uncertainties caused by ripples and bad 
channels. 



^GILDAS is developed and distributed by the Observatoire de Grenoble and IRAM. 
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Results 



3.1. Line identification and measurement 



A total of 74 individual features were detected in our survey, includi ng some less than ce r 



tain detections. Line identification is mainly based on the JPL c atalog (IPickett et al. 



and the Cologne database for molecular spectroscopy (CDMS, iMiiller et al 



2001 



1998 



i 



20051' 



which give the molecular line frequencies based on theoretical calculations. We also used 
the NIST Recommended Rest Frequencies for Observed Interstellar Molecular Microwave 
Transition^, which provides the molecular line frequencies derived from observations of var- 
ious sour c es. M o reover, we utilized the r ecent observations of lRC +10216 and CRL618 by 
He et al.l (120081 ). ICernicharo et al.l (120001 ). and jPardo et al.l (120071 ) for reference. The com- 
plete spectra obtained with the ARO 12 m telescope and the SMT 10 m telescope are plotted 
in Figs. [1] and [2] along with identified lines marked. Figs. [3] and H] give the detected molecular 
line profiles. Note that the two features at 221.459 GHz and 241.458 GHz clearly seen in 
Fig. [2] are the strong CO lines from the image sideband. Among the detected features, 69 
are identified to arise from 21 molecular species and their isotopic variants. Table [1] lists 
their assignments, peak and integrated intensities and line widths species-by species. The 
half maximum line widths [FWHM) were obtained by fitting Gaussian line profiles. If a 
line has an unresolved hyperfine structure, the peak intensity of the strongest transition is 
given, and the integrated intensity and the FWHM are those of the combined feature. 

Five faint features remain unidentified and they are listed in Table [2J Comparing with 
the unidentified line listing in the NIST frequency table, we find that the 148444, 255940, 
and 265936 MHz lines have no previous recorded detections elsewhere, and the 257035 and 
2622 55 MHz may have corr esponding U lines at 257015 and 262273 MHz detected in Sgr B2 
(N) jNummelin et al.lll998h . 



A detailed description of the molecules detection is given below. 

CO. The J = 2-1 transitions of CO, ^^CO, and C^''0 were detected at the 1.3mm 
window. Fig. |4] clearly shows that the shape of the CO line differs from that of the ^^CO 
line. The CO line shows a parabolic profile, while the ^^CO line has a double- peaked profile, 
sugges ting the former is optically thick, while the latter is optically thin. iTeyssier et c 
( 120061 ) observed the CO transitions from J = 1-0 to J = 6-5 in CIT 6. The shape of the 



^http:/ /spec.jpl. nasa.gov. 
http://www.phl. uni-kocln.de/vorhcrsagcn/. 
^http: // physics .nist .gov/cgi-bin/micro /table5 /start .pi 
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CO (2-1) line obtained by our observations is very similar to the ir results (see their Fig. 6) . 
The CO (2-1) and ^^CO (2-1) lines have also been observed by iGroenewegen et al.l (119961 ) 
using the IRAM telescope in 1991. They obtained the integrated intensities of the CO (2-1) 
and ^^CO (2-1) lines to be 230 and 16.5Kkms~^ respectively, a factor of ~ 3 larger than 
our results. As the beam size of the IRAM telescope is about half of that of the SMT, the 
discrepancy of the CO intensity can be explained with different beam dilution effects. Hence, 
we can conclude that the variation of CO intensity is relatively small (within 10%) during 
the past 15 years. The extremely faint C^^O emission should be a new detection for this 
object. 



SiC2- SiC2 has been detected by lWoods et al.l (120031 ) in the spectra of CIT 6 through the 
SiC2 (5o,5-4o,4)) transition at 115.382 GHz. In our hne survey, a number of SiC2 transitions 
were prominently detected at the 2 mm and 1.3 mm windows. These lines have a similar 
profile with an average linewidth of /SVfwhm = 23.9 ± 1.9 km s^^. We also detected two 
faint ^^SiC2 transition although their intensities are only at a 2-3cr noise level. We cannot find 
previous papers reporting on the detections of these SiC2 and ^^SiC2 transitions tabulated 
in Table [H 

CN. CIT 6 is characterized by bright CN emission. The frequencies of three CN (2- 
1) fine-structure groups lie in the region of the SMT spectrum. All of them were clearly 
seen in Fig. ^[ The three CN (2-1 ) fine-structure groups consist of 18 hyperfine structure 
components. iBachiller et al.l (119971 ) carried out a survey of CN (2-1) and (1-0) emission in 
a sample of evolved stars. They detected the strongest two CN (2-1) groups in CIT 6. Their 
observations suggests the intensity ratio of high- and low-frequency fine-structure groups to 
be 1.5, in excellent agreement with our result of 1.6. 

CS. Strong CS (3-2 ] emission was de t ected in th e 2 mm window. Th is line had previ- 
ously been discovered by iBujarrabal et al.l (jl994j ) and iHenkel et al.l (119851 ) . We also clearly 
detected its isotopic transitions, ^^CS (3-2, 5-4) and C^^S (3-2, 5-4). The faint C^^S (5-4) 
emission was only marginally discovered. The CS (5-4) transition at 244.9 GHz does not lie 
in the frequency region observed in the work because of a gap between 244.5-252.5 GHz in 
the survey. Figs. [3] and H] show some evidences that the line profiles of CS and its isotopologue 
emission are slightly asymmetric in shape with a brighter red wing. 

SiS. The J = 6-5 and J = 5-4 transitions of SiS in CIT 6 have previously been 



observed by many researchers (IHenkel et al.lll985l : IBujarrabal et al.lll994j : IWoods et al.ll2003 



Schoier et al.ll2007l ). In this work, three SiS transitions with higher J (8-7, 13-12, and 14-13) 
were clearly detected with an average linewidth of AVfwhm = 22.6 ± 2.2 km s"^. 



SiO. SiO emission in this carbon star has been extensively studied. The transitions 
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observed by previo u s studies include J = 2-1, 3-2, 5-4, a nd 8-7 (iBujarrabal et al.l 11994 
ieging et al.l l2000l : IWoods et al.l l2003l : ISchoier et al.l l2006l ) . In this work, we report the 
detection of the SiO (6-5) transition and its ^^S and isotope lines as well. The SiO (6-5) 
line detected in this survey has a similar width to other SiO transitions detected by previous 
researchers. 



C yH. The N = 1-0 transition of C2H in this object has been detected by lFukasaku et al. 
(I1994J ). The = 3-2 transition lies in our surveyed frequency range. The rotation transition 



is split in six hyperfine-structure lines grouped in three fine-structure groups, all of which 
were detected at the 2 mm window. The two main components (37/2-25/2 and 35/2-23/2) are 
quite strong and have a well-defined profile, which is similar to that of the CS line with 
higher flux in the red wing. The 37/2-25/2 transition is blended with a weak SiC2 line. The 
35/2-25/2 transition is too faint to obtain a reliable intensity. 

HCN. The J = 3-2 transitions of HON and H^^CN were clearly detected. The HCN 
(3-2) transition is the sec ond brightest l i ne aft er the CO (3-2) transition in our survey. This 
line was first observed by iBieging et al.l (I2OOOI ) using the SMT, who obtained an integrated 
intensity of 43. 1 K km s~^ , in excellent agreement with our measurement. Given its line shape, 
the line should be optically thick. There are five favorable vibrationally excited lines of HCN 
present in the surveyed frequency range. Three of them were clearly detected, as illustrated 
in Table [H The weak 1/2 = 2"^^ , 2^^ J = 3-2 transitions are below the 3a noise level. In CIT 6, 
HCN is the only species with detected vibrationally excited transitions. Fig. H] shows that 
these vibrationally excited lines are probably narrower than the HCN (3-2) line. Discarding 
the uncertain detection of the 1^2 = 2^-^, 2^^ J = 3-2 transitions and the z/2 = 1^^ J = 3-2 
transition which is partially blended with the HCN (3-2) line, we obtain the average width 
of the vibrationally excited lines of HCN to be 14.2Kkms~^. The narrow linewidth of the 
vibrationally excited lines suggests that they might arise from the hot inner region with low 
expansion velocity. We do not detect the vibrationally excited lines of H^^CN. 



C3N. iFukasaku et all fll994h and IWoods et all (|2003|) observed the C3N (11-10 a,b) 
transitions in CIT 6. Both lines are relatively weak. Due to fine-structure interactions, 
every rotational transition of C3N is split into doublets of similar intensity. There are 12 
C3N transitions present in our survey range with six at the 2 mm window and six at the 
1.3 mm window. We detected the strongest four C3N transitions (14-13 a,b and 15-14 a,b) 
at the 2 mm window with Tr > 20 mK. The 16-15 a,b transitions at around 158 G Hz have 
a comparable intensity with the other four transitions at the 2 mm window. However, they 
fall within a spectral region with a high noise level, and thus are not listed in Table. [H 



C4H. To our knowle d ge, C 4H has not been detected in CIT 6 before this work. iFukasaku et al. 



( I1994J ) and IWoods et al.l (120031 ) estimated the intensity upper-limits of the C4H (10-9 a,b) 
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transitions, which indicate that C4H emission in this object is relatively weak. Analogous to 
that of C3N, every rotational transition of C4H is split into two components with a similar 
intensity. There are 12 favorable lines of C4H present in the frequency range surveyed here. 
Eight of them were detected with Tji = 9-35 mK. The C4H (15-14 a,b) and (16-15 a,b) 
transitions were clearly detected with well-defined profiles at the 2 mm window. The C4H 
(14-13 a,b) lines fall within a spectral region with a high noise level, and thus are not listed in 
Table [H The four C4H lines (24-23 a,b and 28-27 a,b) at the 1.3 mm window are extremely 
faint, and are only marginally detected. The C4H (25-24 a,b) transitions at around 238 GHz 
are overwhelmed by noise. 

HC3N. Seven HC3N transitions from J = 15-14 to J = 29-28 are in the frequency range 
of our survey. All of them were prominently de tected with Tr > 45mK. Previously only 



lower J transitions of HC.sN have been reported (IHenkel et al.lll985l : iBujarrabal et al. 



1994 



Fukasaku et al.lll994j : IWoods et al.ll2003l ) and to the best of our knowledge all these detections 
are new. Because of a gap in the 1.3 mm spectrum, the HC3N (27-26) at 245.606 GHz was 
not detected here. As shown in Table [H there is a trend that the widths of HC3N lines 
decrease with increasing J-values, suggesting that the high-J transitions might originate in 
hot inner regions which have a lower velocity compared to the regions from where the low- J 
transitions arise. All HC3N lines detected in this survey show asymmetric profiles (see Fig. [3] 
andH]). Opposite to CS and C2H lines, HC3N lines have a brighter blue wing relative to the 
red one. This indicates that the molecular envelope is asymmetric along the line of sight, i.e., 
the chemical compositions and/or physical conditions in the red and blue sides are different. 

CH3CN. There are 14 favorable CH3CN transitions in the survey region. We detected 
the blended features of the (12i-lli) and the (12o-llo) transitions wi t h Tr ~ 23 mK, which 
are the strongest two among the 14 CH3CN transitions. IWoods et al.l (120031 ) failed to detect 
CH3CN in this object and only estimated the intensity upper limit of the (6i-5i) transition. 
This is therefore the first detection of CH3CN in GIT 6. 



All the lines discov ered in this survey have been detected in IRG+10216 (iHe et al.ll2008 



Gernicharo et al.ll2000l ). Gompared to the spectrum of IRG+10216 in the same frequency 
range, the non-detected species in the spectrum of GIT 6 include G3H, G3H2, G2S, G3S, 
H2GO, SiG, SiN, PN, and metal containing molecules, all of which have only weak emission 
in the spectra of IRG+10216, and thus are below our detectior i limit. A number of v ibra- 
tionally excited species have been detected in IRG+10216 (see iGernicharo et al.l l2000l . and 
the references therein). For GIT 6, however, no vibrationally excited lines except those of 
HGN are strong enough to be detected. Moreover, we find no evidence for the presence of 
ionic species in GIT 6. 
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3.2. Rotation diagram analysis and fractional abundances 

The standard "rotation-diagram" method was applied to calculate the excitation tem- 
peratures (Tea;) and column densities (N) of the molecules observed in our spectra. From 
the equation of radiative transfer and assuming that the lines are optically thin, the level 
populations are in local thermal equilibrium (LTE), and T^x » T^g, where T^g is the cosmic 
background radiation temperature (2.7K), we have the well-known relation, 

In — = In — „ ^ = In ■ '^^ 



Nu, gu, and Eu are the population, degeneracy, and excitation energy of the upper level, 
J Tgdv is the integration of the source brightness temperature over the velocity, S is the line 
strength, is the dipole moment, z/ is the line frequency, and Q is the rotational partition 
function. If several transitions arising from levels covering a wide energy range are observed, 
Tex and can be deduced using a straight-line fit to Nu/gu versus Eu/kT^x- 

The rot at ion- diagram provides important tools for studies of excitation conditions. De- 
parture from the linear relation can be caused by different excitation mechanisms or misiden- 
tification. For SiC2, SiS, HC3N, and C4H, there are adequate numbers of detected transitions 
covering a wide range of excitation energy and their rotation diagrams are given in Fig. O 
Good linear correlations were obtained for these species. For the calculations, we have cor- 
rected the effect of beam dilution through = Tr{61 + 9"^)/ 9^, where 9b is the antenna full 
beam at half-power (~ 40" and 30" for the ARO 12 m and the SMT respectively) and 9s 
is the s ource diameter. 9^ ma y be different for different species. We followed the assump- 



tion by iFukasaku et al.l (119941 ) and took an commo n of 20" for all the species. Since the 



source size is likely to vary from species to species (iLindqvist et al.ll2000l ). this assumption 
will introduce a ~ 50% uncertainty in the derived column densities. The derived excitation 
temperatures and column densities are given in Table [31 For these species for which only 
one line was detected or observed transitions arise from a narrow range of energy levels, the 
method of rotation diagram method cannot be employed and a constant T^x of 40 K was 
assumed for the calculations of their column densities. 

Assuming that the molecular envelope is a spherical shell, the emission is optically thin, 
Tex is uniform throughout the envelope, mass loss rate and expansion velocity are constant 
during the formation of the envelope, and the molecular density follows an law, we 
determined the fracti onal abundance s of the observed species with respect to H2 through the 



formula proposed by lOlofssonI (119961 ). 



fx = 1.7 X 10-"^ ^ " ^ , (2) 
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where J T^dv is given in Kkms~^, the full half power beam width 6b is in arc sec, the 
expansion velocity fg is in kms^^, D is the distance in pc, Mn^ mass loss rate in 

Mq yr~^, Uui the line frequency in GHz, gu is the statistical weight of the upper level, A^i is the 
Einstein coefficient for the transition, Ei is the energy of the lower level, and Xj^g = Ri,e/ {(^bD) 
with Ri and Re the inner radius and o uter radius of the shell. For the calculations, we adopted 
D = 400 pc JCohen fc Hitchonlll996h . 



For the determination of molecular abundances using eq. 2, we first derive the expansion 
velocity V(.=18 km s~^ from the profile of the CO (2 -1) line. The rnass lo ss rate of = 
3.2 X 10-6 M© yr-i is obtained by applying eq. 2 of lwinters et all J2002h to the CO (2-1) 
line assuming fco = 1 x 10"^. These values of Ve and M-^^ are in excellent agreement with 



those given by iFukasaku et al.l (119941 ). These parameters suggest that CIT6 has a similar 
expansion velocity and about one order of magnitude lower mass loss rate in comparison 
with IRC+10216. 

The resultant abundances are given in Table. [31 Combined with the uncertainties in- 
troduced by calibration, noise, basehne, and source size, we estimate that the errors of the 
column densities and abundances amount to a factor of ~ 2. One should bear in mind 
that when the emission is optically thick the and fx listed in Table. [3] rep r esent only 
lower limits. A com parison of our results with those derived by iFukasaku et al.l (Il994j ) and 
Woods et al.l (120031 ) are also given in Table. [3] and no significant discrepancies are found. 



Discussion 



4.1. Chemistry 

4.1.1. Oxygen-bearing molecules 



Similar to IRC+10216, the molecular envelope of CIT 6 is also characterized by a lack 
of oxygen-bearing compounds and abundance of carbon-bear ing compounds. In IRC -l-10216, 
three 0-bearing molecules are observed, CO, SiO and HCO+ (jCernicharo et al. I I2OO0I ) whereas 
in CIT 6, only two 0-bearing molecules (CO and SiO) are detected. We find that the SiO 
(6-5)/^^CO (2-1) integrated intensity ratio in CIT 6 is 1.05, in excellent agreement with the 
value of 1.13 in IRC+10216 (iHe et al.ll2008l ). In our discussion, ^^CO is taken as a reference 
molecule because the ^^CO (2-1) line is likely optically thin and the ^^C/^^C isotopic ratios 
are similar in CI T 6 and IRC-t-10216 (see Sect. 4.2). SiO has been commonly detected in 
C-rich envelopes JSchoier et al 1l2006h . implying the presence of icy comets surr ounding the 
stars ( Agundez fc Cernicharol 20061 ) and/or non-equilibrium chemical processes ( Cherchnefi 
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20061 ). HCO"'" has a relatively low abundance in IRC+10216 (iGlassgoldl Il996l ) . and is pre- 
sumably below the current detection limit even if it might be present in CIT6. 



4- 1-2. Carbon-bearing molecules 

We observed abundant carbon chains and radicals in CIT6, including CO, SiC2, CN, 
HCN, CS, C2H, C3N, C4H, HC3N, and CH3CN, all of whi ch are linear. This cha racteristic 
feature is similar to those of IRC+10216 and TMC 1 (see ICernicharo et al.ll2000l ) although 
these lines are much fainter in CIT6. 

The most intriguing characteristic of CIT6 is the strong CN emission. The integrated 
intensity ratio of the CN (2-1) group and the CO (2-1) transition is 4.6, a factor of 2.2 
larger than the value in IRC+10216 (IHe et al.l 120081 ). CN is mainly formed through the 
photodissociation of HCN, 

HCN + /iz/ ^ CN + H. (3) 

According to iHe et all fl2008l ). the H^^CN (3-2)/i3CO (2-1) integrated intensity ratio in 
IRC+10216 is 4.5, a factor of 3.2 larger than that in CIT6. Therefore, our observations 
provide strong evidence that reaction ([3]) dominates the chemistry of CN and HCN in AGB 
stars and the photodissociation is more efficient in the more evolved C-rich envelope CIT6. 

The above discussion also suggests that about 30% CN formed from HCN has been 
destroyed. On the other hand, CN can be reprocessed into HC3N through the reaction 



CN + C2H2 HC3N + H. 



(4) 



We do find that the HC3N line intensities relative to the ^'^CO (2-1) transition in CIT6 are 
a factor of ~ 3 larger than those in IRC+10216, indicating efficient formation of HC3N in 
CIT6. We did not find evidence for the enhancement of the C3N radical, suggesting that 
photodissociation of HC3N into C3N is insignificant in this object. 

CIT6 shows strong C2H emission. The C2H radical is dominantly produced through 
the photodissociation reaction 



C2H2 + hu^ C2H + H. 



(5) 



Our observations show that the C2H line intensities relative to the ^^CO (2-1) transition in 
CIT6 are almost the same at those in IRC+10216, suggesting that there is no significant 
C2H enhancement in CIT6 compared to IRC+10216. This probably has an implication that 
C2H is dominantly processed in early AGB stages. Since C2H and HC3N have the same 
chemical precursor, the /(C2H)//(HC3N) abundance ratio can be used to test the chemical 
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formation path (IWootten et al.lll980l ). Our results yield a ratio of 4.2 for .f(C2H)//(HC.sN) 
which is in good agreement with those found in interstellar clouds by IWootten et al.l (llQSOl ) 
and is consistent with the prediction of gas phase chemistry. 

On the other hand, as shown in Figs. [3] and HI HC3N lines show profiles that differ from 
those of the C2H lines. This suggests that the /(C2H) / /(HC3N) ratio is not a constant in the 
envelope. The chemistry structures in the red and blue sides of the star are not identical. This 
is probably a consequence of an inhomogeneous density distribution or physical environment. 
A complete understanding of the molecular environment of CIT 6 calls for a comprehensive 
3D photochemistry model and high-resolution mapping observations. 



C4H is positively detected in CIT 6. The C4H radical can be formed via 

C2H2 + C2 — > C4H + H. 



(6) 



Nejad fc Millarl (119871 ) suggested that ion-molecule reactions also play an important role 
for the pr oduction of C4H, unlike those for HC3N. Based on millimeter interferometer ob- 



servation, iDayal fc Biegingl (119931 ) found that the photochemical model underestimates the 
C4H abundance in IRC+10216 by a factor of 5. Our observations suggest that CI T 6 has a 



f ('C4 H) //(HC.sN) abundance ratio of 3.1, which is lower than that in IRC+10216 (IHe et al. 



20081 ) by a factor of ~ 3. On the other hand, C4H can be photodissociated into C2H. How- 
ever, as no enhancement of C2H is found in CIT 6, destruction of C4H should be insignificant 
in this C-rich envelope. Therefore, the high radical abundance of C4H in IRC+10216 still 
remains mystery. 

CH3CN is a new finding for this C-rich envelope. This molecule is a symmetric top and 
has been widely used as diagnostics of excitation temperature. We detected only one weak 
CH3CN line in CIT 6, and thus cannot use it to derive an excitation temperature. CH3CN 
can be produced via 



CH+ + HCN ^ CH3CNH^ 



(7) 



followed by 



CH3CNH+ + e ^ CH3CN + H. (8) 

The formation of CH3CN may be very efficient in CIT 6 since strong HCN emission is de- 
tected. CIT 6 has a CH3CN (12-ll)/i3CO (2-1) i ntegrated interi sity ratio of 0.058, in good 
agreement with the value of 0.043 in IRC+10216 JHc et allboosl ). 



Strong CS emission has been detected in our survey. According to lWillacy fc Cherchnefi 



(119981 ) ■ CS can be rapidly destroyed by shocks w hich might occur a fter a star leaves the 
AGB stage and ejects material in a very fast wind (iHerpin et al.ll2002l ). Therefore, the high 
abundance of CS in CIT 6 suggests that shocks are not important for the chemistry in this 
C-rich envelope. 
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4-1.3. Silicon-bearing molecules 



Three refractory Si-bearing species (SiO, SiS, and SiC2) were detected in CIT 6. Al- 
though other Si-bearing species were detected in IRC+10216 (jCernicharo et aLll200Cll ). emis- 
sion from the other Si-bearing species is relatively faint and should be below the detection 
limit of our observations of CIT 6. 

We find th at the abu n dance s of SiO and SiC2 in CIT 6 are similar to those in IRC+10216 
determined by iHe et al.l (120081 ) . Although the situation may be complicated by optical 
depth effects, the HCN/SiO line intensity ratio has the potential to provide a useful tool to 
discriminate between C- rich and 0-rich env elopes and is a good tracer of mass loss rate for 
M and S stars (see e.g. iBiering et al.lboooh . The HCN (3-2) / SiO (6-5) inte nsity ratio in 



CIT 6 is 8.4, in good agreement the value of 9.7 in IRC+10216 (IHe et al.ll2008l ). There is no 
evidence showii ig that the HCN/SiO line intens ity r atio has depen d ance on the mass loss rate 
of C-rich stars. I Gonzalez Delgado et al.l (120031 ) and ISchoier et al.l (120061 ) found a correlation 
between the mass loss rate and the SiO abundance for AGB stars. This is described as 
freeze-out of SiO molecules onto dust grains. The similarity of the SiO abundances in CIT 6 
and IRC+10216 suggests that the depletion of SiO onto dust grains might be insignificant 
for the two C-rich envelopes. 

Our observations show that the SiS abundance in CIT 6 is lower than that in IRC+10216. 
The SiS ( 14-13) /S i O (6- 5) intensity ratio in CI T 6 is 1.1, about half of that in IRC+10216 
found by iHe et al. J2008h . ISchoier et al.l (120071 ) did not find a strong correlation between 
the mass loss rate and the SiS abundance, suggesting that SiS molecules are less likely to 
be depleted onto dust grains than SiO molecules. Therefore, freeze-out should not be the 
reason for the depletion of SiS in CIT 6 since no depletion of SiO is found for this object. On 
the other hand, with a high eff iciency, SiS can be photodis sociated into Si"*", which initiates 
circumstellar SiC„ chemistry (IMacKay fc Charnleyl Il999l ). Hence, we infer that efficient 
photodissociation in CIT 6 has destroyed SiS, and silicon chemistry has been ongoing in 
this evolved C-star envelope. Interferometric observations of these Si-bearing molecules are 
obviously needed to verify the conjecture. 



4.2. Isotopic ratios 

Isotopic ratios of various elements provide substantial tests for nucleosynthesis of low- 
and intermediate- mass stars (LIMS). When a LIMS evolves into the AGB stage, the nucle- 
osynthesized products synthesized through the CNO cycle inside the star are dredged up to 
the surface and then are ejected into the circumstellar envelope. Consequently, the isotopic 
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composition in the circumstellar shell can be markedly changed. Based on the fractional 
abundances proposed in Table [3], we deduce the isotopic ratios (or their lower limits) of 
carbon, oxygen, silicon, and sulfur in CIT6. The results are listed in Table HI The errors 
estimated from the measureme nt and calibratioii are g iven. For com parison, we als o list the 
isotopic ratios for IRC+10216 (ICernicharo et al.ll2000l ) and the Sun (jLodderd 120031 ). 



4-2.1. Carbon 



The ^^C/^^C abundance ratio is the most studied isotopic abundance in LIMS. Stan- 
dard stellar models predict that the ^^C/^^C abundance ratio can be significantly increased 
during the nucleosynthesis and dredge-up processes in the AGB stage. However, extensive 
observations have shown that the ^^C/^^C abundance r atios in LIMSs are considerably lowe r 
than those e xpecte d by standard stellar models (e.g. ICharbonnel fc do Nascimentd Il998l ). 
Charbonnell (119951 ) proposed an extra mixing process to account for the low ^^C/^^C. In 
low-mass AGB stars, the nonstandard mixing called cool bottom processing may dec rease 
the ^^C/^'^C ratio to ~ 4 (jSackmann fc Boothroydlll999l : iBoothroyd &: Sackmannlll999l ). For 
AGB stars more massive than ~ 4:Mq, the hot bottom burning may take place a nd induce 
i2(~;^i3q further decrease to its equilibriur n value of ~ 3.5 (Frost et al.lll998 ). C urrent 



observations of the CO isotopologues in PNs (IBalser et al.l 120021 : iJosselin fc Bachillerll2003l ) 
suggest that the ^^C/^^C ratio is in the range of 2.2-40, supporting the theory including 
nonstandard mixing processes. 

Three ^^C-bearing species have been detected in this survey, including ^^CO, ^^CS, and 
H^'^CN. However, their main lines are likely optically thick. Therefore, the abundance ratios 
of CO, CS, and HCN and their isotopologues only provide l ower limits of the ^ ^C/^^C ratio. 
Our results are in good agreement with those presented by ISopka et al.l (119891 ) who derived 
the abundance ratios i2CO/i3co> 25 ± 10 and Hi2CN/Hi3CN> 5.4 for GIT 6. 

We have detected the rare isotopes, ^^C'^^S and ^^C'^^S. If the '^^S/^^'S abundance ra- 
tio were known, we co uld obtain the ^^C/^'^C ratio using the two optically thin species. 
Cernicharo et al.l (120001 ) found that the sulfur isotopic ratios in IRC+10216 are close to so- 
lar. Therefore, we reasonably assume that the ^^S/^^S ratio in GIT 6 is the solar value. 
It follows that we obtained the ^^C/^^C ratio of 45.4 ± 4.9 in GIT 6, which is in perfect 
agreement with that found in IRG+10216 and is significantly lower than the solar value. 
The ^^G/^'^G ratios found in t hese C-rich e r ivelop es are also lower than the value of 75 ± 9 
in the Orion Bar proposed by iKeene et al.l (119981 ) using the G^'^O/^'^G^^O abundance ratio. 
This is consistent with the hypothesis that nonstandard mixing processes have decreased the 
i2q^i3q ratios in the envelopes around AGB stars. 
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As shown in Table HI different ^^C/^^C values are obtained for CIT6 if different species 
are used for the calculations. If completely ascribing the ^^C/^^C discrepancies found for 
CIT6 to the opacity effects of the main lines, we can estimate the optical depths of the 
CO (2-1), CS (3-2), and HCN (3-2) lines, which are 1.3, 1.2, and 2.0, respectively. If the 
isotopic lines are also optically thick, the optical depths obtained here should be considered 
as lower limits. 



4-2.2. Oxygen 



The nucleosynthesis and dredge-up processes in t he AGB stage can lea d to strong enrich- 
ment of ^^O relative to ^^O and ^^O (see Busso 2006. for a recent review). Wannier fc Sahai 



( 119871 ) found that the ^^0/^*^0 ratios in C-rich envelopes are markedly higher than the ter- 
restrial and interstellar values, but the ^^0/^*^0 ratios are comparable to the solar value. 

We have detected C^^O, allowing us to derive the ^^O/^^O ratio in CIT 6. The C^^O/C^^O 
abundance ratio give a lower limit of 237 ± 26. We may also use the optically thin species 
^^C^'^O and ^^C^'^O to derive the ^"^O/^^O ratio. As s uming ^^c/isq^ 45 4 (^ggg above), we ob- 
tain a i SQy'iTQ j^atio of 890z b97. ICernicharo et al.l (I2OOOI ) did not obtain the oxygen isotopic 
ratios. iKahane et al.l (Il992l ) calculated the ^^O/^^O ratio for a sample of C-rich envelopes. 



They found that ^^0/^'^0= SAO^tn and 840+^30 for CIT 6 and IRC+10216, respectively. 
These values are in good agreement with our result, agree with each other, and are lower 
than the solar value by a factor of about three, which is consistent with predictions of stellar 
models. 

No ^^0-bearing species was detected. The non-detection of the C^^O (2-1) transition 
at 219560 MHz seems to suggest that ^'^0/^^0> 1 in CIT 6. However, this result should be 
taken with some caution since the C^^O line is very close to the edge of the spectrum. 



4-2.3. Silicon and sulfer 



The nucleosynthesis in LIMS is expected to hardly affec t the elements in the 3rd row of 
the periodic table, such as Si and S. ICernicharo et al.l (I2OOOI ) indeed found that the Si and S 
isotopic ratios in IRC+10216 are compatible with the solar values. 

The silicon isotopes, ^^Si and ^°Si, have been detected through faint emission from ^^SiO, 
^°SiO, and ^^SiC2. Since the SiO and SiC2 lines are probably optically thick, the lower limits 
of the ^^Si/^°Si and ^^Si/^^Si ratios were derived. Given the low abundances of ^^Si and 
3°Si, the 29Si/3°Si ratio is not affected by opacity effects. We obtained 29Si/3°Si= 1.0 ± 0.4, 
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comparable with the values in IRC+10216 and in the Sun. 

The C^^S/C^'^S ratio gives a lower limit of ^^S/^^S in CIT6. C'^^S was only marginally 
detected. The two optically thin species, C^^S and C^'^S, give a ^^S/^^S ratio of 0.2 ± 0.2. 
Consequently, we did not find a significant deviation of the S isotopic ratios in CIT 6 from 
those in IRC+10216 and in the Sun. 



4.3. Is IRC+10216 unique? 

Given its brightness and abundant molecular emission, IRC+10216 is the most surveyed 
object and has frequently served as a standard reference for studying circumstellar chemistry. 
A commonly asked question is whether IRC+10216 can truly represent C-rich AGB stars. 
To investigate the problem, we systematically co mpare the spect ra of CIT 6 with those of 



IRC+10216 obtained in the same survey program (IHe et al.ll2008l ). 



We find that all lines detected in this work have also been observed in IRC+10216, and 
all the strong lines detected in IRC+10216 are seen in CIT 6. In Fig. [6l we compare the 
intensity ratios of the lines detected in both objects. The results show that the line intensity 
ratios in the two objects are in agreement within one order of magnitude with an average 
value of 0.08 with a standard deviation of 0.05. The intensity ratios, however, might be 
affected by beam dilution. CIT 6 is a more compact object and thus suffers from a larger 
beam dilution effect. If we correct for the beam dilution effect, the J(CIT 6)//(IRC+10216) 
ratios should increase by a factor of ^, where 

e = (l + ^)/(l + ^^^^)- (9) 
^cne '^IRC+10216 



Assuming ^cit6 = 20" and 6'irc+io216 = 30" (iFukasaku et al.lll994j ) and using 6b = 40" and 



30" for the ARO 12m and the SMT respectively, we obtain the ^ values of 1.8 and 1.6 for the 
12 m and SMT data, respectively. While recognizing that eq. 9 assumes uniform brightness 
temperature which is probably not realistic, the result does suggest that there is no large 
correction factor difference between the A 2 mm and 3 mm bands. Our conclusion that for 
most of the species the /(CIT6)//(IRC+10216) line ratios are in good agreement therefore 
stands. Consequently, we conclude that IRC+10216 is indeed likely representative of C-rich 
envelopes although comparisons with a larger sample of objects will be needed to make a 
stronger statement. The richness of molecular lines in the spectra of IRC+10216 is mainly 
due to its relatively nearby distance, and not due to any special circumstances. 

Fig.[6]also shows that there are a few molecular species for which the /(CIT 6)//(IRC+10216) 
ratios depart from the average value. For the CN and HC3N lines, the ratios are higher. 



-17- 



whereas for HCN, SiS, and C4H, the ratios are lower. As discussed in Sect. I4.1.2[ this 
partly reflects the chemical evolution in the circumstellar envelope around the more evolved 
AGB star CIT6. However, the cause of the abnormally strong C4H emission in IRC+10216 
remains unknown. 

Far-IR spectra of CIT6 and IRC+10216 have been obtained by the ISO Long Wave- 
length Spectrometer (LWS) (jSchoier et al.l 120021 : ICernicharo et al.lll996l ). The rotational 
transitio ns revealed by the ISO spectra can trace the inner regions of the circumstellar en- 
velopes (IHerpin et al.ll2002l ). In Fig [71 we compare the far-IR spectra of the two objects. 
The spectra were retrieved from the ISO archive. Inspection of the figure shows that their 
far-IR spectra are dominate d by thermal co i itinuu m emission from the dust with some su- 
perimposed molecular lines. iLindqvist et al.l (I2OOOI ) fitted the dust continuum of IRC+10216 
with a single blackbody of 510 K and the continuum of CIT 6 with two blackbodies of 1000 
K and 510 K. In the long wavelength region, the dust temperatures of the two objects are 
therefore almost identical. 

Fig. [H] gives the continuum-subtracted ISO LWS spectra of the two C-rich envelopes. A 
nu mber of lines frorn CO, HCN, H^'^CN, and vibrationally excited HCN have been identified 
by lCernicharo et al.l (119961 ) in the spectrum of IRC+10216. Fig.[8]shows that the relative flux 
ratios of the lines detected in the two object are in good agreement although the detection 
of some weak lines is difficult in CIT 6 as it is fainter. The far-IR line fluxes detected in 
IRC+10216 are a factor of about 10 higher than those in CIT 6, consistent with that found 
for the molecular lines at millimeter wavelengths. However, we find that the far-IR dust 
continuum emission in IRC+10216 is a factor of about 20 times higher than that of CIT 6, 
as shown in Fig [3 This probably suggests that IRC+10216 has a larger dust-to-molecular 
gas ratio. 



4.4. Abundance variations expected in later evolutionary stages 



Since CIT 6 has be en proposed to be a h ighly evolved AGB star on the verge to become 
a Proto-PN (or PPN) (ISchmidt et al.ll2002l ). it would be useful to compare the molecular 
abundances of CIT 6 with the corresponding molecular abundances in the evolved AGB star 
IRC+10216 and a PPN. The archetypical PPN CRL618 has abundant molecular emission. 
It is one of the brightest molecu lar sources in the sky and is an ideal target for investigating 
circumstellar chemistry (see e.g. lPardo et al.ll2007l ). All three objects are carbon rich, and are 
commonly assumed to belong to a common evolutionary sequence. Any systematic difference 
in molecular abundances in these three objects can be used to infer on the processes of 
chemical synthesis and destruction. 
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Fig. [9] gives the fractional abundances of the species in the three objects as a function of 
their photodissociation rates taken from the UMIST dat abase for a temperat ure of 300 K. The 



abundances in IRC+10216 and CRL 618 are taken fromlWoods et ah (2003) andlPardo et al. 



(120071 ). respectively. To facilitate the comparison, we follow iPardo et al.l (120071 ) and give the 
molecular abundances relative to HC3N. According to this figure, except for CN and HCN, 
the molecular abundances in CIT6 are obviously closer to those in IRC+10216 than to 
CRL 618. CRL 618 has lower abundances of SiO and CS, suggesting that destruction by 
shocks and depletion onto dust grains may play an important role for the chemistry in this 
object. Compared to the other two objects, CRL 618, having a BO central star, is exposed in 
a stronger UV radiation field. However, we do not find a correlation between the abundance 
differences of these objects and the molecular photodissociation rates, suggesting that the 
destruction by photodissociation is not a factor affecting the chemistry in these objects during 
the AGB-PPN transition. However, we do find that some molecules have been efficiently 
reprocessed during this transition. Since the dynamical timescale of PPN is ~10^ years, our 
results therefore suggest a rapid change of the chemical compositions after the star evolves 
into the PPN stage. 

We should note that this study is limited to simple molecules in the gas phase. There 
is strong evidence from infrared spectroscopy that the solid-state phase chemistry is very 
active in the AGB-PPN evolutionary transition, with many aromatic and aliphatic com- 
pounds being formed in the circumstellar envelope (jKwokil2004l ). Even some simple gas-phase 
molecules (e.g., acetylene and benzene) are difficult to detect via rotational transitions in 
the mm/submm region, and complex organics even more so. Consequently, the technique 
of rotational mm/submm spectroscopy is unable to serve as a complete diagnostic tool for 
circumstellar chemistry. It is more useful as a probe of photo- or shock-chemistry in the PN 
phase, or during the formative stage of early AGB evolution. 



5. Conclusions 

The presence of rich molecular species around evolved stars provides an opportunity to 
study the evolution of chemistry in circumstellar envelopes, which have been widely suggested 
as one of the main sources of organic compounds in space. As part of our project of inves- 
tigating circumstellar chemistry, this paper reports a spectral line survey of the carbon-rich 
envelope CIT6, covering the frequency range between 131-160, 219-244, and 252-268 GHz 
with a high sensitivity. A total of 74 lines are reported in the survey. We identify 69 lines 
belonging to 21 different molecular species and isotopologues, most of which are carbon- 
bearing species. The new species include two carbon-chain molecules, C4H and CH3CN, and 
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seven C, O, S, and Si isotopologues. Several new transitions from known species have been 
detected for the first time in this object. The species with the largest number of detected 
emission lines in our survey is SiC2, which has 19 lines. It is followed by HC3N, with 7 lines. 

We find that the line profiles for some molecular species have different shapes, sug- 
gesting that the chemical structure is asymmetric in the envelope. A comprehensive 3D 
photochemistry model is required to account for the line intensities and profiles in CIT6. 

The excitation temperatures, column densities and abundances of the detected molecules 
are determined through rotation-diagram analysis. The spectra of CIT 6 are characterized 
by a large CN/HCN abundance ratio. Our results suggest that there is evidence for the 
photodissociation of HCN and SiS and the formation of CN and HC3N in the evolved AGB 
envelope. The strong SiO and CS emission may suggest that depletion onto dust grains and 
destruction by shocks are insignificant in this object. An abundance comparison with the 
PPN CRL618 implies to a rapid chemical evolution after a star leaves the AGB stage. 

In order to investigate whether the molecular environment of IRC+ 10216 is intrinsi- 
cally unique, we systematically compare its spectra with those of CIT 6. According to the 
comparison, we find that the molecular species can be classified into three groups, a) for 
most of the species, the intensity ratios of individual lines in the two objects are in good 
agreement with each other; b) the emission from HC3N and CN may be enhanced in CIT 6; 
c) the emission from SiS, HCN, and C4H may be depleted in CIT 6. The differences of the 
line-intensity ratios in the two objects are probably a consequence of chemical evolution with 
the exception of C4H, for which a high abundance in IRC-l-10216 cannot be explained by 
photochemical models. The ISO LWS spectra show that CIT 6 has a lower continuum-to-line 
ratio than IRC-l-10216, suggesting that the latter might have a larger dust-to-molecular gas 
ratio. 

Using the same telescope settings, we also obtained the spectra of the AGB stars 
IRC+10216 and CRL3068, the PPN CRL2688, and the young PN NGC7207. A detailed 
study of the chemical compositions in different evolutionary stages will be published in a 
separate paper. 

We are grateful to the ARO staff for their help during the observing run. Otto Peng 
and Yu-Chin Huang assisted in the observations. We also thank Jun-ichi Nakashima and 
Jin-Hua He for useful discussions. We acknowledge an anonymous referee for comments 
which helped strengthen the paper. The work was supported by a grant from the Research 
Grants Council of the Hong Kong Special Administrative Region, China (Project No. HKU 
7028/07P). 
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Fig. 1. — Spectrum of CIT6 in the frequency range 131-160 GHz obtained with the ARO 
12 m telescope. The spectra have been smoothed to a resolution of 1 MHz. 
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Fig. 2.— Spectrum of CIT 6 in the frequency ranges 219-244 GHz and 252-268 GHz obtained 
with the SMT 10 m telescope. The spectra have been smoothed to a resolution of 3 MHz. 
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Fig. 2. — continued. 
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Fig. 3. — Line profiles of the molecular species detected in the ARO 12 m spectrum. The 
spectral resolution is 1 MHz. The labeled hues are shown at zero velocity in each panel. 
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Fig. 4. — Line profiles of tlie molecular species (ietected in the SMT spectrum. The spectral 
resolution is 3 MHz. The labeled hues are shown at zero velocity in each panel. 
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Fig. 4. — continued. 
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Fig. 5. — Rotational diagrams for the detected species in CIT6. Error bars are given on the 
lower right. 
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Fig. 6. — Integrated intensity ratios of the lines detected in CIT6 and those detected in 
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Fig. 7. — Upper panel: the ISO LWS spectra of IRC+10216 and GIT 6, for which the fluxes 
were normalized such that F^Qnm — 1; Lower panel: the absolute- flux ratio of the two spectra. 
Note that the ratio at A > 160 //m bears large errors. 




Fig. 8. — The continuum-subtracted ISO LWS spe ctra of IRC+10216 and CIT6. The sohd 
hues mark the CO and HCN emission identified by lCernicharo et al.l (Il996l ) in the spectrum 
of IRC+10216. 
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Table 1. Molecular transitions detected in CIT6. 



Frequency 


Species 


Transition 


rms 




jTndv 


AVpwHM 


Remarks" 


(MHz) 




(upper-lower) 


(mK) 


(K) 


(K km/s) 


(km/s) 




o o r' o o r\ 

230538. U 


CO 


J=2— 1 


4.2 


2.555 


67.590 


23.7 




220398.7 




J=2-l 


5.6 


0.163 


4.996 


29.0 




224714.4 


C 'O 


J=2-l 


7.1 


0.022 


0.258 


25.6 




137180.8 


SiC2 


6o,6-5o,5 


6.8 


0.040 


1.040 


24.2 




140920.1 


SiC2 


62,5-52,4 


4.8 


0.062 


1.459 


22.6 




141751.5 


SiC2 


64,3-53,2 


6.1 


0.073 


1.671 


23.4 




141755.4 


SiC2 


64,2-54,1 












145325.8 


SiC2 


62,4-52,3 


6.5 


0.088 


1.867 


21.1 




158499.2 


SiC2 


7o,7-6o,6 


9.5 


0.086 


2.548 


25.3 




220773.7 


SiC2 


10o,io-9o,9 


5.6 


0.087 


2.282 


22.8 




222009.4 


SiC2 


92,7-82,6 


6.4 


0.069 


1.887 


24.8 




232534.1 


SiC2 


102.9-92.8 


4.0 


0.072 


1.937 


25.2 




234534.0 


SiC2 


108,2-98,1 


4.8 


0.030 


0.452 


16.1: 




234534.0 


SiC2 


108,3-98,2 












235713.0 


SiC2 


106,5-96,4 


5.2 


0.045 


1.101 


25.2 




235713.1 


blC2 


106,4-96,3 












237150.0 


DIO2 


104,7-94,6 


5.8 


0.038 


0.859 


20.7 




2i7i6l.6 


&1O2 


104^6-94,5 


5.8 


0.036 


0.771 


23.5 




241307. 7 


bl02 


llo,ii-10o,io 


4.2 


0.074 


1.770 


22.7 






aiC2 


112,10-102,9 


3.9 


0.094 


2.448 


23.6 




zooUbo.U 


bl02 


ll8,3-10s,2 


9.3 


0.039 


0.471 


10.5: 




2580d5.0 


SiC2 


118,4-108,3 












259433.3 


bl02 


116,6-106,5 


7.8 


0.054 


1.301 


27.5 




259433.3 


S1O2 


116,5-106,4 












2d1150.7 


blC2 


114,8-104,7 


4.1 


0.044 


1.049 


21.5 




261509. 3 


blC2 


114.7-104,6 


4.1 


0.051 


1.270 


26.5 




2d199U.7 


SiC2 


12o,12-ll0,ll 


7.1 


0.387 


9.945 


26.1 
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29o:r^ 
"blC2 


llo,ii-10o,io 


5.3 


0.013 


0.366 


35.3 




zOoDoi.O 




11.,., 1 n., _ 
Il6,6-J^06,5 


A 


n ni n 




o2.y. 


* 


255631.8 


29SiC2 


116,5-106,4 










* 


226287.4 


CN 


N(J,F)=2(3/2,l/2)-l(3/2,l/2) 


4.9 


0.034 


1.278 


53.0: 


* 


226298.9 


CN 


N(J,F)=2(3/2,l/2)-l(3/2,,3/2) 










* 


226303.0 


CN 


N(,T,F)=2(3/2,3/2)-l(3/2,l/2) 










* 


226314.5 


CN 


N(,T,F)=2(3/2,3/2)-l(3/2,3/2) 










* 


226332.5 


CN 


N(,T,F)=2(3/2,3/2)-l(3/2,5/2) 


4.9 


0.052 


1.819 


25.6 


* 


226341.9 


CN 


N(J,F)=2(3/2,5/2)-l(3/2,3/2) 










* 


226359.9 


CN 


N(J,F)=2(3/2,5/2)-l(3/2,5/2) 










* 


226616.6 


CN 


N(J,F)=2(3/2,l/2)-l(l/2,3/2) 


4.9 


0.053 


1.593 


26.8 


* 


226632.2 


CN 


N(J,F)=2(3/2,3/2)-l(l/2,3/2) 










* 


226659.6 


CN 


N(J,F)=2(3/2,5/2)-l(l/2,3/2) 


4.9 


0.178 


6.092 


29.3 


* 


226663.7 


CN 


N(J,F)=2(3/2,l/2)-l(l/2,l/2) 










* 


226679.3 


CN 


N(J,F)=2(3/2,3/2)-l(l/2,l/2) 










* 


226874.2 


CN 


N(J,F)=2(5/2,5/2)-l(3/2,3/2) 


4.9 


0.394 


12.176 


24.6 


* 


226874.8 


CN 


N(J,F)=2(5/2,7/2)-l(3/2,5/2) 










* 
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Table 1 — Continued 



Frequency Species Transition rms Tr / T^dv AVpwHM Remarks" 

(MHz) (upper-lower) (mK) (K) (K km/s) (km/s) 



226875.9 


CN 


N(J,F)=2(5/2,3/2)-l(3/2,l/2) 


— 


— 


— 


— 


226887.4 


CN 


N(J,F)=2(5/2,3/2)-l(3/2,3/2) 


— 


— 


— 


— 


226892.1 


CN 


N(J,F)=2(5/2,5/2)-l(3/2,5/2) 


— 


— 


— 


— 


226905.4 


CN 


N(J,F)=2(5/2,3/2)-l(3/2,5/2) 


— 


— 


— 


— 


146969.0 


cs 


J=3-2 


4.8 


0.303 


7.058 


21.6 


138739.3 


13CS 


J=3-2 


4.6 


0.023 


0.495 


21.0 


231221.0 


13CS 


J=5-4 


5.0 


0.027 


0.846 


27.6 


242913.6 


c-«s 


J=5-4 


7.4 


0.032 


0.563: 


12.4: 


144617.1 


c^^^s 


J=3-2 


4.9 


0.040 


1.039 


26.1 


241016.1 


C34S 


J=5-4 


4.6 


0.062 


1.315 


23.1 


145227.0 


SiS 


J=8-7 


6.5 


0.067 


1.745 


22.1 


235961.1 


SiS 


J= 13-12 


5.2 


0.119 


2.598 


20.6 


254102.9 


SiS 


J=14-13 


3.6 


0.213 


5.791 


25.0 


260.518.0 


SiO 


J=6-5 


4.8 


0.233 


5.233 


20.3 


254216.7 


3«SiO 


J=6-5 


3.6 


0.022 


0.581 


21.5 


257255.2 


2«SiO 


J=6-5 


6.9 


0.021 


0.583 


30.3 


262005.3 


C2H 


Nj=37/2-26/2 


7.1 


0.387 


9.945 


26.1 


262066.1 


C2H 


Nj=35/2-23/2 


7.1 


0.244 


6.596 


25.4 


262208.4 


C2H 


N,T=35/2-26/2 


6.0 


0.024 


0.5.39: 


22.9: 


265886.4 


HCN 


J=3-2 


3.7 


1.959 


43.942 


21.2 


265852.8 


HCN 


1/2=1'^" J=3-2 


3.7 


0.067 


1.392 


20.4 


267109.1 


HCN 


i/2=22/ J=3-2 


4.2 


0.011 


0.098: 


22.5: 


267120.1 


HCN 


i/2=22« J=3-2 


— 


— 


— 


— 


267199.3 


HCN 


V2=l^f 3=3-2 


4.2 


0.050 


0.876 


14.4 


267243.2 


HCN 


1/2=2° J=3-2 


4.2 


0.020 


0.318 


14.0 


259011.8 


Hi^CN 


J=3-2 


9.6 


0.296 


7.120 


20.8 


138515.7 


C3N 


N=14-13 a 


7.3 


0.039 


0.747 


21.5 


138534.5 


C3N 


N=14-13 b 


7.3 


0.031 


0.674 


25.5 


148409.1 


C3N 


N=15-14 a 


4.2 


0.021 


0.537 


29.9 


148427.8 


C3N 


N=15-14 b 


4.2 


0.021 


0.521 


35.6 


142728.8 


C4H 


N=15-14 a 


6.0 


0.022 


0.567 


33.7 


142767.3 


C4H 


N=15-14 b 


6.0 


0.021 


0.450 


23.1 


152243.6 


C4H 


N=16-15 a 


5.7 


0.018 


0.433 


23.6 


152282.1 


C4H 


N=16-15 b 


5.7 


0.020 


0.505 


27.2 


228348.6 


C4H 


N=24-23 a 


4.1 


0.015 


0.224: 


19.7: 


228387.0 


C4H 


N=24-23 b 


4.1 


0.035 


0.256: 


13.1: 


266389.9 


C4H 


N=28-27 a 


3.9 


0.008 


0.140: 


29.9: 


266428.2 


C4H 


N=28-27 b 


3.9 


0.009 


0.175: 


19.1: 


136464.4 


HC3N 


J=15-14 


6.6 


0.160 


4.131 


24.0 


145560.9 


HC3N 


J=16-15 


6.5 


0.139 


3.772 


24.1 


154657.3 


HC3N 


J=17-16 


7.0 


0.131 


3.165 


24.4 


227418.9 


HC3N 


J=25-24 


6.3 


0.064 


1.404 


22.0 


236512.8 


HC3N 


J=26-25 


6.2 


0.053 


1.078 


21.7 


254699.5 


HC3N 


J=28-27 


3.9 


0.068 


1.160 


21.7 
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Table 1 — Continued 



Frequency 


Species 


Transition 


rras 






AVp-wHM 


Remarks'" 


(MHz) 




(upper-lower) 


(mK) 


(K) 


(K km/s) 


(km/s) 




263792.3 


HC3N 


J=29-28 


7.9 


0.048 


0.990 


19.6 




220743.0 


CH3CN 


12(1)-11(1) 


5.6 


0.023 


0.290 


17.6: 


* 


220747.2 


CH3CN 


12(0)-11(0) 










* 



"♦-unsolved hyperfine structure lines; $-blended with other species. 



Table 2. Unidentified lines. 



Frequency 


rms 


Tr 


/ T^dv 


^Vfwhm 


(MHz) 


(mK) 


(K) 


(K km/s) 


(km/s) 


148444 


4.5 


0.008 


0.249 


14.8 


257035 


7.1 


0.015 


0.442 


38.6 


255940 


4.9 


0.012 


0.207 


29.4 


262255 


6.0 


0.022 


0.501 


26.7 


265936 


3.7 


0.012 


0.214 


37.3 
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Table 3. Excitation temperatures, column densities and abundances with respect to H2". 



Species 


TexCK)^ 


N{ 


jm ^) 


fx 




This paper" 


F94'' 


This paper'^ 


W03<^ 




1 .0 






A( R\ 


6.l(-b) 


29c:p 




Of^/' 1 
Z.Z0[^LO ) 




3.6(-7) 




olo 


zy.D 






3-4(-bj 


3.1(-bJ 


^4x1 


00. y 


Q HA ^^ A\ 


< 0.04(^io ) 


4.U(^-D ) 


< 1.7(-bj 




41. U 






1.3(-b) 


i.4(-b) 






i.Do(^i ( ) 












1.35(16) 












6.89(14) 








cs 




1.72(14) 




2.0(-6) 


2.5(-6) 






1.28(13) 




1.5(-7) 








6.18(12): 




7.2(-8): 








2.58(13) 




3.0 (-7) 




SiO 




2.15(13) 




7.0(-7) 


1.0(-6) 


30SiO 




2.45(12) 




8.0(-8) 




29SiO 




2.43(12) 




7.9(-8) 




HCN 




1.10(14) 


5.13el5 


1.2(-6) 


l.l(-5) 






1.83(13) 




2.0(-7) 


3.0(-7) 


CN 




1.22(15) 




2.6(-5) 


2.0(-5) 


C2H 




7.43(14) 


4.0 ±2.0(14) 


5.4(-6) 


6.9(-6) 


C3N 




2.16(13) 




1.9(-6) 


2.6(-6) 


CH3CN 




7.88(12) 




l.l(-7) 


< 1.3(-7) 



"'x{y) represents x x 10^; 

''A constant excitation temperature of 40 K was assumed for the species for which 
the rotation diagrams cannot be obtained; 

''For the species with optically thick emission (e.g. CO and HCN), this gives the 
lower hmits; 

'^F94: from Fukasaku et al. (1994); W03: from Woods et al. (2003). 
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Table 4. Isotopic abundance ratios. 



Isotopic ratio 


CIT6 
Sepcies 


Value*^ 


IRC+10216'^ 


Solar^ 




12034g/13032g 


45.4±4.9'^ 


45±3 


89 




i2C0/"C0 


12.1±1.3" 










13.4±1.7" 








Hi^CN/Hi^CN 


6.0±0.7'^ 






160/170 


130160/120170 


890±97-^ 




2680 




C^^o/Ci^o 


237 ±26" 






29Si/30Si 


29SiO/=^0SiO 


1.0±0.4 


1.45±0.13 


1.52 


28Si/=^°Si 


28SiO/30SiO 


8.8±1.9" 


20.3±2.0" 


29.9 


28Si/29Si 


28SiO/29SiO 


8.9±2.8" 


15.4±2.0" 


19.6 




28SiC2/2^SiC2 


6.7±3.1^ 






32g/34g 




6.7±1.2'= 


21.8±2.6 


22.5 


33S/34S 




0.2±0.2 


0.18±0.1 


0.18 



"Prom Cernicharo et al. (2000); 
*Prom Lodders (2003); 

"^The errors are estimated based on the measurement and calibration 
uncertainties; 

"^Assume that the ^^S/^^S ratio is solar; 

"Should be treated as lower limits due to opacity effect; 

^Adopted: i2c/i3c=45.4. 



